Phonon Softening and Direct to Indirect Bandgap Crossover in Strained Single Layer 

MoSe 2 
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Motivated by recent experimental observations of Tongay et al. [Tongay et al, Nano Letters, 
12(11), 5576 (2012)] we show how the electronic properties and Raman characteristics of single 
layer MoSe2 are affected by elastic biaxial strain. We found that with increasing strain: (1) the E' 
and E" Raman peaks (Ei g and E2 g in bulk) exhibit significant red shifts (up to ~30 cm -1 ), (2) 
the position of the A\ peak remains at 180 cm -1 (Ai g in bulk) and does not change considerably 
with further strain, (3) the dispersion of low energy flexural phonons crosses over from quadratic to 
linear and (4) the electronic band structure undergoes a direct to indirect bandgap crossover under 
~3% biaxial tensile strain. Thus the application of strain appears to be a promising approach for a 
rapid and reversible tuning of the electronic, vibrational and optical properties of single layer MoSe2 
and similar MX2 dichalcogenides. 

PACS numbers: 81.16.Pr, 68.65.Pq, 66.30.Pa, 81.05.uc 
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I. INTRODUCTION 

The discovery of graphenei 2 - is expected to play an 
important role in future nanoscience and nanotechnol- 
ogy applications. Recent advances in nanoscale growth 
and mechanical exfoliation techniques have led not only 
to the fabrication of high-quality graphenes but also the 
emergence of several new classes of two-dimensional (2D) 
structures such as ultra-thin transition metal dichalco- 
genides (TMDs). Just like graphene TMDs have hexag- 
onal crystal structure composed of layers of metal atoms 
(M) sandwiched between layers of chalcogen atoms (X) 
with stoichiometry MX2. TMDs have electronic semi- 
conducting or metallic properties and most of bulk 
TMDs possess indirect band gaps of the order 1-2 
eV £ The synthesis of several TMDs has been realized 
experimentally^— and the stability and electronic prop- 
erties of various single layer dichalcogenides was re- 
cently reported^ Recent studies of several semiconduct- 
ing TMDs such as MoS 2 , WS 2 , MoTe 2 and MoSe 2 have 
shown that the band gap increases and transforms to a 
direct band gap with decreasing number of layers^iSr— 
These sizeable band gaps make them well suited for elec- 
tronic applications such as transistors, photodctectors 
and electroluminescent devices. Another unique feature 
of two-dimensional ultra-thin materials is the possibil- 
ity to apply large reversible elastic strain. It has been 
shown that graphene can be strained up to 20 % of its 
ideal structure with only small changes in its electronic 
band structure, which is in contrast to TMDsji 2 - - — 

The most recent efforts have been directed towards the 
synthesis and manipulation of molybdenum disclcnide 
(MoSe 2 ) single layers. In addition to preliminary reports 
on synthesis of few-layer MX 2 structures , 16 ' 1 Tongay et 
al. demonstrated that a single layer of MoSe 2 possesses 
a direct optical gap of 1.55 eV and exhibits good thermal 



stability^ However, to our knowledge, no research ex- 
ists addressing the question of how electronic properties 
and lattice dynamics of single layer MoSe 2 are affected 
by strain. In the present work, we investigate the effect 
of biaxial strain on the electronic and vibrational prop- 
erties of single layer MoSe 2 structures. Our calculations 
revealed that applying biaxial strain is able to tune the 
Raman characteristics and electronic band structure of 
single layer MoSe 2 . We expect that this study will offer 
new opportunities for strain-engineered nanoscale opto- 
electronic device applications. 



II. COMPUTATIONAL METHODOLOGY 

Calculations of physical properties of equilibrium and 
strained structures were carried out in the framework 
of density functional theory (DFT), using the plane- 
wave self-consistent field (PWSCF) code as implemented 
in the QUANTUM-ESPRESSO packaged The gener- 
alized gradient approximation (GGA) of Pcrdcw-Burkc- 
Ernzerhof (PBE) was used for the exchange-correlation 
potential^ A plane-wave kinetic energy cutoff of 40 
Ry and density cutoff of 400 Ry were used. Bril- 
louin zone integration was performed using a shifted 
27x27x1 Monkhorst-Pack meshi^I To eliminate the in- 
teraction emerging from periodic boundary conditions 
calculations are performed with a large unitcell includ- 
ing 12 A vacuum space between adjacent MoSe 2 single 
layers. Ground states and total energies of all systems 
were obtained after full geometry relaxation with forces 
on the atoms smaller than 0.02 eV/A. Phonon frequen- 
cies and phonon eigenvectors were calculated in a 4x4x 1 
q-grid using the density functional perturbation theory 
(DFPT)^ 2 . 



2 




r m k r 




m r KM 



FIG. 1: (Color online) (a) Phonon dispersion for single layer 
MoSe2 and displacement directions of Raman-active modes. 
Inset: top view of atomic structure of MoSe2 (b) Electronic 
band structure and partial density of states calculated within 
GGA and GqWq. The energies are relative to the Fermi level 
(i.e. Ef = 0) Inset: Comparison of the band edges calculated 
by GGA, GGA+SO and LDA. 



III. RESULTS AND DISCUSSION 

A. Single Layer MoSe2 

Atomically thin M0SC2 is the most recently synthe- 
sized member of the ultra-thin transition metal dichalco- 
genides. A single layer structure of MoSe2 can be viewed 
as a 3-layer stacking of Mo and Se atoms wherein molyb- 
denum atoms are sandwiched between layers of trigonally 
arranged selenium atoms. In this configuration, known 
as 1H structure, each Mo atom is coordinated to six Se 
atoms. In the equilibrium geometry of the single layer 
MoSe2, the Mo-Se distance, the Se-Se distance and the 
lattice constant of the hexagonal unit cell are calculated 
to be 2.528, 3.293 and |a|=3.321 A, respectively. Thus, 
in comparison with single layer M0S2 having lattice con- 
stant 3.12 A, MoSe2 has a larger atomic structure due to 
the larger atomic radius of the Se atoms. In Fig. []Ja) 
we present the calculated phonon dispersion and vibra- 
tional characteristics of R-active modes of single layer 
MoSe2. It is seen that the frequencies of all phonon 
branches in the whole Brillouin Zone have positive values 
i.e., MoSe2 crystal can remain stable by generating the 
required restoring force against atomic distortions. Al- 
though the synthesis of single layers of MoSe2 have been 
achieved only on various substrates so far, energy opti- 



mization and phonon dispersion calculations show that 
freestanding MoSe2 monolayers are quite stable. 

The calculated band structure shown in Fig. QJb) 
shows that the single layer MoSe2 is a semiconductor 
with a direct bandgap at the K high symmetry point. 
The top of the valence band is mainly composed of Mo- 
d( x 2_ y 2-) orbitals and, albeit small, also by M.o-d( xytyZtZX ) 
and Se-p( Xjy ) orbitals. Though the electrons occupying 
Ad and 5s shells of a neutral Mo atom are treated as va- 
lence electrons, 5s states have a negligible contribution 
even away from the Fermi level. However, it is seen that 
all Mo-d, Se-p( Xj y) and Se-p( 2 ) orbitals are hybridized to 
form the conduction band edge. Note that x and y com- 
ponent of the orbitals have the same energy due to the 
two-dimensional lattice symmetry. Due to the presence of 
another minimum in the conduction band edge between 
the r and K points we also used local density approx- 
imation (LDA) for the exchange correlation functional. 
Tongay et aZ.— demonstrated that such a neighboring 
conduction band dip is quite sensitive to temperature- 
induced strain and layer-layer interaction. The lattice 
constant of a single-layer MoSe2 is found to be smaller 
(3.22 A) in LDA because of the overestimation of the 
strength of the covalent bonds. The inset of Fig. \V[b) 
shows a zoom of both LDA, GGA and GGA+SO band 
structures along the T-K-M path in the Brillouin zone. 
We see a splitting of the band edges at the vicinity of 
K point when including spin-orbit (SO) interactions. In- 
terestingly, LDA predicts that single layer MoSe2 is an 
indirect semiconductor where the top of the valence band 
is at the K point, and the bottom of the conduction 
band is found to be between the T and K points. It 
appears that the physically correct description of single 
layer MoSe2 structure that was reported experimentally 
as a semiconductor— with a direct bandgap is achieved 
by using the GGA functionals. Therefore, in the rest of 
our study PBE calculations will be employed. 

Although our PBE calculation captures the qualita- 
tive nature of the band edge states it underestimates 
the band gap. To correct the band gap we use the 
Go Wo approximatio n 24-26 as implemented in the VASP 
package i 27 i 28 It was shown before that this kind of cal- 
culation has to be converged with respect to the vacuum 
spacing! 29 ! 30 We calculate the quasi-particle shifts using 
(12x12x1) k-point grid with vacuum spacings of 15 A, 
20 A, 25 A, and 30 A and then extrapolate them to infin- 
ity in a similar fashion as done in Ref. (2j|. Our Go Wo 
calculation results in a band gap value of 2.33 eV. Note 
that, excitonic effects are pronounced in 2D materials as 
M0SC2 and we need to calculate the binding energy of the 
exciton in order to estimate the optical gap^2r— Due to 
the parabolic shape of the conduction and valence bands 
we assume that the exciton has the Mott-Wannier char- 
acter and the effective mass approximation can be used. 
Hole and electron effective masses calculated from the 
band structure arc m* h = 0.65 and m* = 0.53 (in units 
of electron mass). Moreover, we account for the 2D po- 
larizability of the system using the model introduced by 
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FIG. 2: (Color online) (a) Phonon dispersions for 1% and 7% 
strained single layer MoSe2 (b) Evolution of the Raman-active 
modes as a function of applied strain. 



Cudazzo et al£^ Using parameters reported by Ref. [3l[ 
we find the binding energy of the exciton 0.7 eV, which 
is in agreement with the value obtained by solving the 
Bethe-Salpeter equation . 31 ' 33 Combining this with our 
result from Go Wo calculation we estimate the optical gap 
to be 1.63 eV, which is close to the reported experimental 
value of 1.56 eV. 



B. Lattice Dynamics and Phonon Softening 

While bulk MoSc 2 crystal has the D 6 /j point group 
symmetry, it turns into D3^ for single layer structure. 
Just like in MoS2r^ lattice vibrations of single layer 
MoSe2 is characterized by nine phonon branches includ- 
ing three acoustic and six optical branches. Analysis of 
lattice dynamics shows that the decomposition of the vi- 
bration representation at the T point is V = A[ + AE 1 + 
IE" + 2A' 2 '. While LA and TA acoustic branches have 
linear dispersion, the frequency of the out of plane flexu- 
ral (ZA) mode has a quadratic dispersion in the vicinity 
of q=0. Near to the T point the in-plane sound velocity 
of the LA and TA modes are found to be 1.5xl0 3 and 
0.9 xlO 3 m/s, respectively. It is also seen that differing 
from graphene's phonon dispersion the acoustic and op- 
tical modes are well-separated from each other. Differing 
from M0S2 and WS2 structures in which the A[ mode is 



the highest Raman-active mode, in MoSe2 the A[ mode 
is located in between the E' and E" modes. In Fig. 
[TJa) the atomic displacements of the R-active modes are 
presented. Interestingly, although recent Raman spec- 
troscopy measurements of Tongay et al*£ revealed that 
only one characteristic Raman peak of single-layer M0SC2 
Eig is observable at 240 cm -1 , our symmetry analysis 
predicts the existence of two more Raman-active modes. 
Significantly decreased intensity of Raman peak around 
290 cm -1 can be explained by considering the coupling 
of in-plane counter-phase vibration of MoSe2 with the 
Si02/Si substrate. Moreover, in experiment the disap- 
pearance of E" (E\ g in bulk) Raman shift around 170 
cm -1 is due to the inactivity of this vibrational mode to 
normally incident light. 

Recent experiments on graphene and other few-layer 
materials have revealed that the application of biaxial 
and uniaxial tensile stress is possible by using flexible 
substrates. We will investigate now how the lattice dy- 
namics of single layer MoSe2 is affected by biaxial strain. 

Fig. [5] shows phonon dispersion curves and the evo- 
lution of Raman peaks as a function of applied strain. 
Below 1% strain, the E' branch of single layer M0SC2 
experiences ~10 cm -1 softening which is larger than 
that of the E g mode (diatomic corresponding to E') 
of graphenei^r— It is seen that the low-energy flexu- 
ral phonons (ZA) of unstrained single layer MoSe2 turns 
almost in a linear dependence under 7% strain. Such 
a change in phonon dispersion implies a significant de- 
crease in scattering of electrons by the flcxural phonons 
and an increase in electron mobility for strained MoSe2. 
For the E' mode there is a small splitting at the T point, 
as in slightly polar materials such as M0S2 and WS2»^ 
It is also seen that due to the weakening of the covalent 
bonds under strain the E' and E" modes (that corre- 
spond to the Ei g and E 2g in bulk MoSe2, respectively) 
exhibit significant red shifts of about ~30 cm -1 . Espe- 
cially, at the K point, the highest optical mode dips even 
below the band corresponding to the E' mode. Another 
remarkable point is the A! 1 mode (corresponding to the 
A\ g in bulk) that vibrates perpendicular to the applied 
strain which experiences a smaller softening and is al- 
most pinned to 231 cm^ 1 for strain values higher than 
1%. Similar direction-dependent softening behavior was 
observed for various types of materials! 36 ' 37 



Strain-Induced Direct-to-indirect Bandgap 
Crossover 



In addition to strain-induced changes in the lat- 
tice dynamics, the electronic properties of single layer 
MoSe2 may exhibit surprising deviations. Existence of a 
bandgap even in the monolayer structure shows strong 
bonding nature and the covalent character of MoSe 2 lay- 
ers. Though the covalent bonds formed from valance 
electrons of Mo and Se atoms are not broken by elas- 
tic strain, the strength of the bond is weaken with the 
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FIG. 3: (Color online) Evolution of the electronic band dis- 
persion of single layer MoSe2 as a function of strain. Fermi 
level shown by red line is set to zero in the energy spectrum. 
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increasing distance. Therefore band energies and disper- 
sion can be expected to be strongly dependent on tensile 
strain. Evolution of the electronic band dispersion un- 
der finite biaxal strain (e = 1 — 7%) is shown in Fig. 
[31 It is clearly seen that the band edges at the T and 
K symmetry points are influenced significantly by biax- 
ial tensile strain. Up to 3% strain, single layer MoSc 2 
remains a direct bandgap semiconductor with both the 
top of the valence band and the bottom of the conduc- 
tion band located at the K point. However, when ^3% 
strain is applied, as a result of the band shifting at the 
K point, two valence band edges with the same energy 
appear in the Brillouin Zone while the conduction band 
edge remains fixed in its original position. Here, with 
~3% strain one has both a direct (K— >K) and an indi- 
rect (r — !>K) bandgap. 

Such an interesting electronic structure may result in 
the presence of two types of holes having different ef- 
fective masses and the coexistence of direct and indirect 
band transitions in the optical spectra. We also found 
that increasing further the strain results in the separa- 
tion (and shifting) of the uppermost valence band to- 
wards the conduction band accompanied by the lowering 
of the conduction band dip at the K point. Moreover, the 
second conduction band edge located (between T and K) 
very close to the dip of the conduction band at the K 
point moves upwards in energy. In addition, band edges 
of M0S2 also experiences similar effects with increasing 
tensile strain. 12 i 13 The evolution of direct (K->K) and 
indirect (K— > TK and T — >K) energy bandgaps under 
tensile stress are presented in Fig. HJa). Both direct and 
indirect bandgaps decrease monotonically with increas- 
ing strain, however the rate of decrease is faster for the 
indirect bandgap. Note that, the variation of band gaps 
with increasing strain has the same qualitative behavior 



FIG. 4: (Color online) Evolution of (a) direct and indirect 
energy band gaps, (b) electron and hole effective mass, nor- 
malized to unstrained value (ml/ml and m£/m£ ), under 
biaxial tensile strain, (c) Energy bandgaps calculated by us- 
ing PBE+SO and PBE+GoWo are shown by solid and dashed 
lines, respectively. Red, green and blue colors are used for 
transitions K— >K, K— >■ TK and V — >K, respectively. 



for PBE, PBE+SO and PBE+G W calculations. We 
predict that the direct-to-indircct crossover takes place 
at 3-4% strain. 

As a consequence of the strain-dependent transition of 
the valence band edge from the K to T point the effective 
mass of the holes (and electrons) propagating through the 
single layer MoSe2 lattice changes dramatically. From 
Fig. [5J it is seen that while the mass of the electrons at 
the band edge does not change notably, the hole mass 
shows significant deviations under strain. Since the va- 
lence band edge is sharpened by biaxial strain the effec- 
tive mass of the holes is decreased and after the bandgap 
crossover holes belong to the uppermost valence band at 
the r point inversely affected by increasing strain. Note 
that the holes under 7% strain are ~2 times heavier than 
that for unstrained single layer MoSe2. 

In summary, motivated by the recent experimental 
study of Tongay et al*& on single layer MoSe2 we in- 
vestigated the electronic properties and lattice dynamics 
of single layer MoSe 2 as a function of biaxial strain. We 
showed that GGA describes more accurately the elec- 
tronic structure of single layer MoSe2 and that LDA 
predicts the wrong nature of the bandgap. We found 
that monolayer MoSe2 has four R-active characteristic 
modes: among these, the modes having counter-phase 
in-plane motion arc significantly modified by strain and 
they soften by increasing strain. Red-shift in the position 
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of these resonant Raman peaks are predicted to be ~30 
cm -1 . However the A[ branch (corresponds to the A lg in 
bulk) is negligibly influenced by strain. Moreover the lin- 
ear dispersion of flexural mode in strained-MoSe2 implies 
a reduced scattering of electrons from these ZA phonons. 
Our findings also revealed that the electronic band struc- 
ture of single layer MoSe2 also undergoes substantial 
changes under biaxial tensile strain. While the direct 
bandgap linearly decreases up to ~3% tensile strain, a 
direct to indirect bandgap crossover occurs when strain is 
further increased. Moreover the appearance of a positive 
curvature at the K point in the phonon spectra and pos- 
sible metallicity for higher strain values suggests super- 
conductivity in highly strained single layer MoSe2. Our 
findings indicate that the bandgap of single layer MoSe2 



can be tuned reversibly by biaxial strain and is able to 
capture a broad range of the solar spectrum. 
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